Bound states of interacting helium atoms 
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We study the possibility of existence of bound states for two interacting 4 He atoms. It is 
shown that for some potentials, there exist not only discrete levels but also bands akin to 
those in the Kronig-Penney model. 

1. INTRODUCTION 

It is known that most of the paradoxes of modern microscopic theory of superfluidity [1] 
have to do with the commonly accepted assumption of the single-particle Bose condensate 
(SBC) prevailing in the quantum structure of the superfluid component. In the 4 He Bose 
liquid, because of strong interaction, the SBC is strongly suppressed (it comprises no 
more than 10% of all 4 He atoms [2]), consequently it cannot be the basis of the superfluid 
component. If in some domain of momentum space there is strong enough attraction 
between atoms, bound pairs of bosons can form, whence a pair coherent condensate 
appears, and it is this condensate that will be the basis of the superfluid component. 

The hypothesis of coupling of helium atoms below the A point may resolve some of 
those paradoxes [1,3]. 

The main purpose of this paper is to explore the possibility of existence of bound states 
for two interacting 4 He atoms for different theoretical and phenomenological interatomic 
potentials. We will show that for some of the potentials, discrete levels do exist. Besides, 
we will analyze certain "quasi-crystallic" models, motivated by the conjecture that at 
temperatures below the A point, helium possesses a structure close to that of quantum 
crystals. Specifically, we will consider a model of atoms "smeared" over a sphere and a 
Kronig-Penney-type model. These models also allow for pairing of helium atoms; more- 
over, the Kronig-Penney-type model leads to a band of allowed energies. 

2. TWO INTERACTING ATOMS 

During several decades, different phenomenological potentials of interatomic inter- 
action in gaseous and liquid 4 He have been constructed, based on empirical data on 
thermodynamic, kinetic and quantum mechanical properties of helium [4-8]. 

Consider the Schrodinger equation for two helium atoms (\1/ being the radial part of 
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the wave function, L = 0) 

$(r) being the potential and m the reduced mass (2m = mn e = 6.6466 • 10~ 24 g). For a 
given $(r), we analyze Eq. ([]]) for discrete levels by means of solving it by the Runge- 
Kutta method from the left and from the right and trying to choose e such that the 
logarithmic derivatives of the two solutions match at some midpoint. If this fails, it is 
concluded that there is no level. 

Table 1 

Energy levels for two helium atoms assuming various interaction potentials 



Potential 


$(r), 10 12 erg (r in A) 


Level, K 


Rosen-Margenau [5] 


925exp( 4.4r) 560exp( 5.33r) ±4^ _3 


None 


Slatter-Kirkwood [6] 


770exp( 4.6r) ±4^ 


None 


Intem-Schneider [8] 


1200 exp( 4.72r) ±# ±# 


None 


Lennard- Jones 


Case 1 [7]: a = 2.556 A, e = 10.22 K 
Case 2 [9]: a = 2.642 A, e = 10.80 K 


None 
0.0201 


Buckingham [10] 


1 


f 770exp( 4.6r) ^ , r < 2.61; 
977exp( 4.6r) ^ 2^1 ; r > 2.61 


0.00632 


Massey-Buckingham [10] 


1000exp(-4.6r) - ^ 


0.0622 


Buckingham-Hamilton [10] 


977exp(-4.6r)-i^- 2 ^i 


0.00229 



3. AN ATOM AT THE CENTER OF A LATTICE 



At temperatures below the A point, the structure and dynamical properties of helium 
are close to those of quantum crystals (e.g., Putterman [11] pointed out a possibility of 
transverse oscillations in superffuid helium). Motivated by this, let us study the behavior 
of a helium atom in the center of a lattice. For qualitative purposes, we replace the actual 
cubic lattice, with the lattice constant d = 4.5 A, by c = 8 atoms "smeared" over a sphere 
of radius a = d a/3/2. The potential corresponding to such a model is 



U(r) = - $(Vr 2 + a 2 - 2ra cos 6) sin 9 dO , 
2 Jo 



(2) 



where $ is taken to be the Lennard- Jones potential (case 2 above). We then have to find 
a level in the potential U(r), depicted below. 

Obviously, since U (r) is infinite at r = a, there will be infinitely many levels. We find 
the ground level to be E = —11.45 K, and the (certainly unphysical) first excited level to 
be E = +35.64 K. 




Fig. 1. The ground level for the atom in the center of a lattice. 
4. TWO NEIGHBORING CELLS AND A BAND 



Now, consider a "quasi-crystallic" model. A more realistic potential, without an infi- 
nite barrier, is 



!$(r), 0<r<a-d, 
$(a — d), a-d<r<a + d, 

$(2a-r), a + d<r<2a. 



(3) 



The parameter d is taken to be 2.2 A, which corresponds to twice the quantum chemical 
radius of the helium atom. 

We get two levels: an even one with E = —13.862 K and an odd one with E = 
-13.8486 K. 
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Fig. 2. The potential and the levels (the splitting is not noticeable) in the model of two 

neighboring cells. 



Finally, consider a periodic potential: 
$ p (r) = $ 2 (r), 0<r<2a; $ P (r + 2a) = $ P (r) . 



(4) 



We find the allowed band by a well-known method [12]: The wave function is represented 

as 



ip(r) 



AKr 



u K (r) , 



(5) 



where uk{t) is a periodic function. Matching the logarithmic derivative leads to 



cos 2Ka = 



[m(0K(q) + m(oK(o)] - MQK(q) + u 2 (a)u' 1 (o)] 

2[«i(0)«' 2 (0) - « 2 (0)«i(0)] 



(6) 



where tti(r) and ^(r) are two arbitrary linearly independent solutions of the Schrodinger 
equation with energy E. The values of E for which this equation can be solved for K (that 
is, for which the right-hand side is not greater than 1 by absolute value), form an allowed 
band. Solving the Schrodinger equation numerically and plotting the E dependence of 
the right-hand side, one finds the band and its boundaries at different values of d. The 
maximal width of the band, at d ~ 2.34, is about 1 K. Depicted below is the potential 
and the band at d — 2.33. 
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Fig. 3. The allowed band in the periodic potential model. 
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When d is further increased, the upper edge of the band will "overflow" the potential 
and merge into the continuous spectrum. In fact, the band edges coincide with the even 
and odd levels in the two-well problem considered above. The reason is that the even 
and odd (with respect to a) wave functions of that problem (the boundary conditions for 
which are if)(0) = if) (2a) = 0) can be extended into the antiperiodic (cos2fTa = —1) and 
periodic (cos2i^a = 1) wave functions, respectively, of the periodic potential. According 
to the above, these functions correspond to the band edges. 

For a rectangular potential (Kronig-Penney model [13]) of wells of width a divided by 
barriers of height Vq and width b with the lattice constant being fixed: a + b = 4.5 A, the 
band is determined by the equation 

/3 2 -a 2 

— - — — sinh f3b sin aa + cosh f3b cos aa = cos k(a + b) , (7) 
zap 

where a = \ / 2mE/h, f3 = ^j2m(V — E)fi. The values of E for which the right-hand side 
is between -1 and 1, form the allowed band. 



Table 2 

Results (E\ ow , £high and AE — lower edge, upper edge, and band width, respectively, 
all energies are in K): 



a 




E\ ow 




AE 


2.5 


10 


2.15 


2.35 


0.20 


3.5 


10 


1.17 


1.62 


0.45 


3.5 


5 


0.78 


1.56 


0.78 


3.5 


2 


0.38 


1.52 


1.14 


4.0 


2 


0.20 


1.92 


1.72 



5. CONCLUSION 

Thus, we have shown that in a superfluid Bose liquid, formation of bound pairs of 
bosons is possible - - for some potentials, there are discrete levels. Moreover, it has 
been shown that bound states can also exist in more realistic models — the one of atoms 
"smeared" over the sphere and the periodic Kronig-Penney type model. Note in particular 
that the Kronig-Penney type model allows for the existence of bands. The presence of such 
bands may explain the observable critical velocities and the energy gap for the excitation 
spectrum curve ("S curve") in superfluid helium, which will be discussed in subsequent 
papers. 

We are grateful to P.I. Fomin for suggesting the problem and constant encouragement 
and to A. Kostyuk for an interesting discussion and advice. Calculations have been made 
using Mathematica. 
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